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Abstract 

We study resonances of multidimensional chaotic map 
dynamics. We use the calculus of variations to deter- 
mine the additive forcing function that induces the largest 
response, that is, the greatest deviation from the unper- 
turbed dynamics. We include the additional constraint 
that only select degrees of freedom be forced, correspond- 
ing to a very general class of problems in which not all of 
the degrees of freedom in an experimental system are ac- 
cessible to forcing. We find that certain Lagrange multi- 
pliers take on a fundamental physical role as the efficiency 
of the forcing function and the effective forcing experi- 
enced by the degrees of freedom which are not forced 
directly. Furthermore, we find that the product of the 
displacement of nearby trajectories and the effective total 
forcing function is a conserved quantity. We demonstrate 
the efficacy of this methodology with several examples. 



1 Introduction 

Sinusoidally driven nonlinear oscillators have been 
widely studied in contexts ranging from synchroniza- 
tion (l3|) and nonlinear response phenomena S B to 
stochastic resonance (0; (H) and nonlinear transport phe- 



nomena ( 14). Somewhat less studied are system iden- 
tification via resonance curves of nonlinear systems (Q) 
and periodically driven chaotic systems dl lb . An area that 
has received much less attention is resonance phenomena 
of nonlinear systems due to aperiodic and chaotic forc- 
ing functions (8). Plapp and Hiibler (10) and others hi 
have used the calculus of variations to show that a spe- 
cial class of aperiodic driving forces can achieve a large 
energy transfer to a nonlinear oscillator Such nonsinu- 
soidal resonant forcing functions yield a high signal-to- 
noise ratio which can be used for high-resolution system 
identification In a recent paper, Foster, Hiibler, and 
Dahmen lil) explored resonant forcing of chaotic map dy- 
namics in which every degree of freedom in a multidi- 
mensional system is forced. 

In this paper, we present a methodology for determin- 
ing the resonant forcing of a multidimensional chaotic 
map in which only select degrees of freedom are forced. 
This is motivated by the difficulty or impossibility of forc- 
ing all of the degrees of freedom in certain experiments. 
For example, consider the forced one dimensional Henon 
map with a delay: x^"+^^ = l-a{x^"^f +cbx^"-^^ +F''"\ 
This can be written as a two dimensional system with no 
delay, but then only one of the two dimensions is forced, 
corresponding to the presence of an F^"' term but the 
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absence of an F^"^'' term. Therefore, the method we 
present may be appHed to a very general class of prob- 
lems. In practice, the calculations necessary to determine 
the optimal forcing are simpler than in the case where all 
degrees of freedom are forced, particularly if only one de- 
gree of freedom in a high dimension system is forced. We 
show analytically that the resonant forcing functions are 
closely related to the unperturbed dynamics of the system 
in that the product of the displacement of nearby trajecto- 
ries and the effective total forcing function is a conserved 
quantity. Furthermore, we find that certain Lagrange mul- 
tipliers take on a fundamental physical role as the effi- 
ciency of the forcing function and the effective forcing ex- 
perienced by the degrees of freedom which are not forced 
directly. We demonstrate the efficacy of the methodology 
with several examples. 

2 Resonant forcing of select degrees 
of freedom 

We begin with the iterated map dynamics 

x("+i) =f(xW)+FW, (1) 

where x'") G R'^ denotes the state of the t/-dimensional 
system at the «th time step, with n ^Q,l,...,N — 1 and 
f(") g (jgjjotes the forcing function at time step n. This 
system has d degrees of freedom. We define the total forc- 
ing effort to be the magnitude of F: 



N-l 



(2) 



Given the corresponding unperturbed system y("+^) = 
f(^y(")j with y^"' = x'''^ we define the final response as 
the deviation from the unperturbed dynamics 



(3) 



We require that < < d degrees of freedom be un- 
forced. Without loss of generality, we choose to order the 
variables so that xi , . . . .x^^^ are unforced and x^/^^^i , . . . .xj 



are forced. Thus we will require that 

for / = 1,. . . ,du and n = 0,1, 



,N-1. 



(4) 



The Lagrange function L used to determine the forcing 
function that produces the largest response R is 



L = 




(5) 



where A , { 7*"' , . . . , 7^,"' } , and { , , . ■ ■ , M^"^ } are 
Lagrange multipliers and F is a constant. We seek sta- 

in) 

tionary points of L corresponding to dL/dx\ = and 

in] 

dL/dF-' = for all n and / = 1,. . . ,d. These equations 
of motion yield: 

af(") + - 



= 
- 0, 



(6) 
(7) 



where J^j 



{dfi/dxj 



is the Jacobi matrix evaluated 



atx'"'. We have also defined the vector r'"' ^ifj'Lx 7/"'^j' 
where is the unit basis vector in the direction of x,. The 
superscript T indicates the transpose operator. For x'^^^ 
we have the additional equation 



xW-yW+M'"-''=0. 
We now define the quantity 

q(") = p(") +r'"'; 



(8) 



(9) 



after we eliminate the vector Lagrange multipliers 
{/i'"), . . . j/i'^^'' }, the equations corresponding to the 
stationary points of Eq. (|5]l reduce to the same form as 
when all variables are forced [See (0)]: 



g(«) 

ag(^-') 



(10) 

(11) 



Accordingly, we identify G as the effective total forcing 
function; it reduces to the optimal forcing F when we re- 
move the constraint in Eq. (|4]l. We identify the Lagrange 
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multipliers {7, , . . . jy^ } to be the effective forcing expe- 
rienced by the degrees of freedom j for which — 0; 
this changes the trajectories of these degrees of freedom 
via the coupling in f rather than direct additive forcing 
via F. The control is stable if, on average, the displace- 
ment of nearby trajectories decreases. Consider a tra- 
jectory given by Eq. ([T]i, and a nearby trajectory given 
by x'"+'' = f(x'"'') +F''"\ where x and x are related by 
g(n) ^ x'"' — x'"'. If we Taylor expand f(x) for small e, 
we obtain 

£<"+'' =j("'eW. (12) 

Multiplying both sides of the transpose of Eq. ( fTol ) by e , 
we have 

(G("+i))V+''e("+'' = (GW)^e("+i). Using Eq. 
this becomes 

(G("+>ye("+2) = (GW)^e("+>), (13) 

a quantity that is invariant for all n. We define this to be 
the conserved quantity P: 

G(o).e(i)=p = GW.e(«+i)^ (14) 

and note that P depends on the observables x and F as 
well as the Lagrange multipliers F, which we have iden- 
tified as the effective indirect forcing of certain degrees 
of freedom. This further reinforces the idea that G repre- 
sents the effective forcing experienced by the system, tak- 
ing into account the coupling via f. As in the case where 
all degrees of freedom are forced, P is conserved even if 
the unperturbed dynamics is chaotic or periodic. 

2.1 Resonant forcing functions with small 
magnitude 

For weak forcing, we can iterate Eq. ^ and Taylor ex- 
pand for small F. We obtain [see Eq. ( |56] | in the Ap- 
pendix] 

xW-yW^MG^'^-^^-il, (15) 

where we have defined M = / + 
L^r/ 7(^-1) • • • • • (7(^^-1))^ where 



7 is the identity matrix and 

^ ^ pfW-l) ^ j(A'-l)p(A'-2) _^ . . . ^ . • • J('))f(''^ 

(16) 

Using Eq. ( fTTT l. we obtain 

MG("-')-n-AG('^-'). (17) 

Eqs. (|2|i, (|4|i, ([TT] ). and (fTOl i form a complete system 
of equations whereby the unknown forcing amplitudes 
{pf^^ , fI'^ } with i ~ {dii + 1), ■ ■ ■ ,d and Lagrange 
multipliers A and | yj"' , . . . , yj^ } with j = 1, . . . ,du can 
be uniquely determined. At this point it is possible to 
write the final response in terms of the total effort: 

R^ = XF^- (18) 

the details of the proof are in the Appendix [see Eq. (|58])1. 
Thus A is the effective efficiency /F^ of the forcing 
function. One approach to solving this system of equa- 
tions is to treat Eq. (fTTT l as an inhomogeneous eigenvalue 
problem. To find the optimal forcing, we first find the 
eigenvectors {vi,...,Vc/} and eigenvalues A, of the ho- 
mogeneous problem (M — A7)G''^^'' = 0. Then we can 
build the inverse Hx = [M — Xiy^ : 

d y .y^ 

^a = Lt^, (19) 

in which A is not yet determined. Then the solutions to 
the inhomogeneous problem are 

G^f^-^) ^HxQ.. (20) 

Eq. ([Tol l can then be used to build a set of equations 

GW ^ . . . n = 0, . . . ,A?-2 

(21) 

whereby A, F^"^ and F^"^ can be determined for all n with- 
out matrix inversion. This method allows a natural con- 
nection to the case in which all degrees of freedom are 
forced, that is, du ^ in Eq. (|4]l. For this case, Q. = 
and Eq. ( fTTI ) reduces to the homogeneous eigenvalue 
problem (M- A7)F('^-i) = with solutions F'""') = 
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{vi , . . . , each corresponding to an eigenvalue A,. For we can write the solutions to Eq. ( fTTI l: 
the homogeneous case, if we set A = max {A,} to be the 
largest eigenvalue of M then Eq. ( fTSl l holds and the La- 
grange multiplier A assumes the same meaning as in the 
inhomogeneous case, namely, the effective efficiency of 
F. 



X = -{\ 



JO) 



'1 



3 Examples 

3.1 Resonances of coupled shift maps 

We consider the mapping function for coupled shift maps: 




mod(flx|"'' -l-fciCj"'') 
mod(ax2"'' -t-fccj"') 



(22) 

(n) 

and require that only F2 be forced, that is, — for 
all n. Accordingly, F^"' ~ /'"^Ci, where we have defined 
7 = 7i . Since the Jacobi matrix 



a k 
k a 



(23) 



is symmetric and constant, the eigenvectors of M for any 
are v± = {^^j^Y We denote the two eigenvalues of 
M by A±, with A+ > A_. Because of the nature of the 
eigenvectors, we can always write M explicitly: 



M = 



X- 

A- 



(24) 



where we have defined the quantities A* = (A+ ± A_) /2. 
It can be shown that for — 2'' with b e Z+, the eigen- 
values of M ai-e A± = nil \l + {a±k) \. As an ex- 
ample, we set = 2 so that the eigenvalues of M are 
A± = 1 + {a±kf. Likewise, we find i2 = 
using Eq. ( fTSI ). From Eq. ^ we impose the additional 
constraints Fj'*^' = Fj''' =0. If we define 



/3± = ±a/(1 +fl2)2 + 2/t2(fl2 - 1) +;^, 



(25) 



r 



1) 



1 

2k' 
1 

2ak 
1 



3k^ + P±)F. 



'1 



1) 
2 ' 

1) 
2 ' 



= — -(l-a^-k^- 
2a^ 



(1) 



(26) 
(27) 
(28) 
(29) 



Using the normalization condition in Eq. (|2]l, we can de- 
termine the magnitude of /v 



(1) 



^2 — 



2aF 



4a^ + (l-a^--k2 + l5±) 



(30) 



We use Eqs. ( fTsT i and ( |26] | to find the final response 



F2 



(31) 



From this equation it can be shown that j3+ gives the min- 
imum response while j3- gives the maximum response. 
For the special case — > 0, corresponding to the case 
of two uncoupled shift maps where one is forced but 
the other is not, the response reduces to the simple form 
R^/F^^ \+a^. We may also compare the effectiveness 
of only forcing X2 to the effectiveness of forcing both x\ 
and X2- This corresponds to removing the constraint in 
Eq. (|4|l and having yj"' = for all /' and n in Eq. (|5]l. 
We will mark solutions of this system with a (~) to avoid 
confusion with the selectively forced system. Thus G = F 
and Eq. ^ reduces to MF*^-') ^ AF*'^-'). We solve 
this system to find A as well as the individual components 
of F for each timestep with N = 2: 



A = l + {a±k), 
F['^^-{a±kfF^'\ 

f['^^±f^'\ 



Fi'^^{a±kYF^\ 



^2 — 



2 + 2{a±k) 



(32) 
(33) 
(34) 
(35) 
(36) 
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Note that Fj ^ 7'"' ; the effective indirect forcing of x\ 
when only X2 is forced is not simply equal to the optimal 
direct forcing Qix\ when both x\ and x^ are forced. Using 
Eq. (fTSb . we can calculate the final response: 

^=A = l + (fl±^)', (37) 

where as before, the (+) sign in front of k corresponds 
to a maximum final response for a > 0, while a (— ) in 
the same position corresponds to a minimum final re- 
sponse. We will assume a > and use (+) henceforth for 
these results. This quantity also reduces to 1 when 
k ^ 0. To compare the effectiveness of forcing only one 
degree of freedom to forcing both, we assume that the to- 
tal forcing effort F is the same in both cases, then use 
Eqs. (ISTT l and (l37T i to obtain the ratio of final responses; 

^,_R^ _l _ 2[l + {a + kf] 

^ - A l+fl2 + fc2 + ^_- (^'') 

Note that I sls k ^ 0. We plot Eq. ^ in Fig. [l] 

for k = 0.3000. Notice that > 1 for all a; that is, for 
the same total effort F, there will always be a greater fi- 
nal response if all the degrees of freedom are forced. This 
is reasonable, given that when both xi and X2 are forced 
there is a more uniform flow of energy into the system 
than when only X2 is forced and xi is passive. However, 

oc 1 for small a and ^ 1 for large a. Since R^ is 
on the same order of magnitude as R^, for certain exper- 
imental situations it may be sufficient (and presumably 
easier) to build an apparatus in which only one degree of 
freedom is forced rather than all. For a < 0, we choose 
the (— ) sign before k in Eq. (|37] |. The result is the same; 

> 1 for small negative a but E^ ^ 1 for a — > —00. 
Furthermore, for this system we can explicitly verify 
Eq.d; (G(i))^e(2) = {G'^°^fe^'^=P. Using Eq. (O, 
we can write this as 

(G(i))^j(i)e(i)^(G(o))^e(i). (39) 

Using Eqs. (|28j and for the left hand side of Eq. ([39]) 



1.5 



0.5 




Figure 1 ; The ratio of final responses E^ versus the param- 
eter a for the same total effort F = 0.001 and k = 0.3000. 
For a > and k ^^0, the final response is greater when 
both degrees of freedom are forced (corresponding to 
E^ > 1). Note that E^ ^ 1 as a ^ 00 or ^ 0. 

we obtain 

17(1) 

(G(i))^j(i)e(i) = + «2 _ 3^2 + (1) 

^ ' 2ak 

+ k{l-a^-k^ + l5±)e^^^] =P. 

(40) 

Using Eqs. ( l27b and ( |29l ), for the right hand side we find 

(G(0))^e(i) = zIlI +«2 _ 3^2 ^ (1) 
^ ' 2ak ^ ' 

+ k{\-a^-k^ + ^±)e^^^]=P. (41) 

Thus for the coupled shift maps with N = 2 we are able 
to analytically verify that P is a conserved quantity. Now 
consider a coupled shift map with a = driven by a forc- 
ing function that is described by Eqs. ( |29] l and (l30l l. The 
response curve as a function of a can by found by starting 
with Eq. (O and iterating Eq. ([ij; 

R^ = k^[F^y+[a,F^'^+F^Y^ (42) 

where Fj"' ™d Fj'^ are functions of a as given by 
Eqs. (|29] l and (|30] |. respectively. We plot Eq. (l42T l as a 
function of a in Fig. |2] In the same figure we plot the re- 
sults of a numerical simulation in which the final response 
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Figure 2: Shift map resonance curve for final response 
/F^ versus the parameter a. The total effort is F = 
0.001 and the other control parameters are qq = 0.6000 
and k = 0.3000. The maximum of the curve is at a = 
qq within machine precision. The solid line indicates the 
analytical expression given in Eq. (|42] | and the triangles 
indicate the results of a numerical simulation with x'f ^ — 
jCj"^ = 0.1000 as the initial condition. 



is found by using Fj"^ and as the forcing functions 
and iterating the system for two time steps. The initial 
condition used is x 
is F = 0.001. 



(0) „ ^(0) 
1 



0.1000 and the total effort 



3.2 Resonances of the one dimensional 
Henon map with delay 

The forced Henon map with delay = 1 — a(x("')^ + 

cfo("^'' _l_ ^(n) written as the equivalent two- 

dimensional system 



bxf 



1 



(43) 

in) 

and require that only be forced, that is, — for 
all n. Accordingly, F*^"' ~ /'"'ei, where we have defined 
7 = 7i . The Jacobi matrix is 



h 



(44) 



while for N = 2 the matrix M is given by 



M 



-2abx. 



(1) 



-2abx^2^ l+c2 + 4fl2(x^'')^ 

JO) .^0)12 



(45) 



where ~ ^2 1 + '^-^^i [■"^2 ] ■ Using Eq. (fTSI l. 
we find £2 = 



' ^ in-,)- We now define 



a = 2a 



l+cxf^~a{xfy 



j3± = ±Jb^ + 2b^{a^ - 1) + (1 + a2 



(46) 
(47) 



As with the coupled shift maps, solving Eqs. (fTOl i and 
([TtI i. and (HI simultaneously yields the following: 



X = -{l+b^ + a^ + j5±), 



JO) 



AO) 



(1) 



cF 



2ba 
1 



'1-b' 

u2 



a 

v2 



(1) 



= -il-b^-a^-P,)F, 



(48) 
(49) 
(50) 

(51) 



Using the normalization condition in Eq. (|2]l, we can de- 
termine the magnitude of F2^^: 

^(1) 2aF 



4a2 + (l-/72_a2_^ 



(52) 



We can also use Eqs. (fTSl l and ( |48] l to find the final re- 
sponse 



R 



1 



a 



(53) 



From this equation it can be shown that j3+ gives the maxi- 
mum response while j3_ gives the minimum response pro- 
vided fl > 0. Now consider, as with the coupled shift map, 
a Henon map described by Eq. (|43j, only with a oq. 



Then we can find the final response as a function of 

4 



(0) 



and ^2 by starting with Eq. ^ and iterating Eq. ([T]l: 



7;2=/,2[^(0)]2^|^(l)_^^^(0) 



r(0) 



-2 + 2cx^P^ — 2a 



where Fj"^ and Fj'' are functions of a as given by 
Eqs. (ISTT l and (|52T i. respectively. We plot Eq. (|54] l as 



(54) 
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a function of a in Fig. [3] Since we use the approx- 



imation ~ ^2'' ™ '^l- *S2}^ we are able to solve 
Eqs. ( [Tol l. ([TTT i. and (|4]i analytically. With this approx- 
imation and with aq = 1.1000, the has a maximum 
at a = 1.0991. If instead we substitute the exact ex- 



pression ' = 1 — ^[xj*''']^ + cxj"' + F^^ into Eq. (l45T l. 
it is possible to solve Eqs. (fTOl i. ([TtT i. and (|4|i numeri- 
cally for ^2"'. Then we recover the expected maximum 
at a = 1.1000. Furthermore, we find that for F = 0.001, 
at maximum the numerical solution of the exact system 
gives /F^ =6.8215, while the analytical approximate 
solution gives R^/F'^ — 6.8118. This is consistent with 
the claim that the methodology presented in this work will 
give the maximum response for a fixed effort. Since any 
forcing other than the exact solution should be subopti- 
mal, we expect a smaller final response when any approx- 
imations are used. 



4 Conclusion 




6.84 



6,82 



6.80 



6.78 



□ V 



0.9 



.0 1.1 1.2 1.3 
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We study resonances of forced multidimensional chaotic 
map dynamics. We constrain the total forcing effort to be 
fixed [see Eq. (|2|i] and seek the forcing function which 
produces the largest response [Eq. (O], subject to the 
additional constraint that certain degrees of freedom are 
not directly forced [Eq. ©I. To determine this forc- 
ing function, we seek the stationary points of the La- 
grange function [Eq. (|5]l] and thereby obtain equations 
which determine the dynamics of the forcing function [see 
Eqs. ([Tol l and (fTTTil. From these equations we identify the 
effective total forcing to be a vector comprising the direct 
forcing and the Lagrange multipliers that represent the ef- 
fective indirect forcing of certain degrees of freedom [see 
Eq. (|9|l]. We demonstrate that the product of the effec- 
tive forcing and the displacement of nearby trajectories is 
a conserved quantity [Eq. (fT4lll. In the case of small forc- 
ing, we show that another Lagrange multiplier represents 
the efficiency of the forcing function [see Eq. (fTSl)]. In 
the limit that we set the number of unforced degrees of 
freedom to be zero, all of the results reduce to the homo- 
geneous case. The methodology presented can be applied 



Figure 3: Henon map resonance curve showing final re- 
sponse R^ jF'^ versus the parameter a. The total effort is 
F — 0.001, the other control parameters are aq = 1.1000 
and k — 0.3000, and the initial conditions are Xj*^' ~ x"^^ — 
0.1000. The solid line indicates the analytical expres- 
sion given in Eq. ( l54l i and the triangles indicate the re- 
sults of a numerical simulation of the simplified system 
in which the Jacobi matrix of the driven system is ap- 
proximated by the Jacobi matrix of the unperturbed sys- 
tem. The boxes indicate a numerical simulation of the 
exact system. The approximate analytical result and the 
approximate numerical simulation curves have a maxi- 
mum at fl = 1.0991. The exact numerical simulation has 
a maximum at a = 1 .1000, which is the expected peak for 
this system, (a) shows the full range of the experiment, 
< a < 6.0, while (b) shows in detail the region near the 
maximum for all three curves, (b) shows that the exact 
result has a greater maximum response than the approxi- 
mate results, as expected. 



7 



to a very general class of problems in which not all of the 
degrees of freedom in an experimental system are acces- 
sible to forcing. 

We demonstrate the effectiveness of the methodology 
with several examples. We compare forcing one degree of 
freedom in a system of two coupled shift maps to forcing 
both degrees and show that the final response is greater 
but on the same order of magnitude when both are forced 
[see Fig.[r|. We present a resonance curve for the coupled 
shift map in Fig. |2] and verify explicitly that the optimal 
effective forcing complements the separation of nearby 
trajectories [see Eqs. (|40] | and WH]. We also apply this 
method to a forced one dimensional Henon map with a 
delay [Eq. (|43])1. a problem which cannot be solved us- 
ing the homogeneous case in which all degrees of free- 
dom are forced. We solve for the optimal forcing function 
analytically for two time steps by approximating the Ja- 
cobi matrix of the forced system by the Jacobi matrix of 
the corresponding unperturbed system. We also solve this 
system numerically without approximations and demon- 
strate that the exact solution reproduces the correct peak 
in the resonance curve and has a greater final response at 
maximum [see Fig. O. Thus we show that the method 
may be used for system identification. In the future we 
plan to compare the effectiveness of this methodology for 
system identification to that of other methods such as pe- 
riodic driving and coupling a test system to a virtual 
model with tunable parameters (j5|). 
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6 Appendix: Proofs 



Proof of Eq. (O: Using Eq. ^ and 



p(A'-n) ^ ^j(A'-n+l)jr|^p(A'-n+l) ^plW-n+ljj _ p(A'-n) 

^ (^j{N-n+l)y I |^j(A'-n+2)-J 7" |-p(A'-«+2) _|_ p(A'-n+2)j p(A'-n+l) | 

+ (j(^-"+l)j^p(A'-n+l) _p(A'-n) 
+ (^j(^-"+l)~)^(^j(W-«+2)~)^p(A'-n+2) 

_ ^j(A'-A'+l)-J^p(Ar-«+l) _|_ ^jlW-n+O-J^plAf-n+l) _ plW-n) 
= (^j(^-n+l)^^(^j(W-n+2)^^ |^p(A'-n+2) _)_p(A'-n+2)j _ plW-n) 



p(A'-«) ^ (•j(iV-n+l))^...(j(A'-l)^^|-p(A'-l)_,_p(A'-l)] _ p(A'-«) 



Using this result. 



_y{N) ^-p{N-l) _^ jiN- 1 )p(iV-2) _j_ j(iV- 1) j(A'-2)p(A'-3) 

+ ...+j(^-i)...j{i)f«') 

+ • • • jC) (jC))^ • • • (j(^-i))^ [f(^-i) +r(^-i)] 

_ |p(iV-l) _|_ j(iV-l)p(A'-2) + . . . + [j(^-l) , 



r(i) 



]p{0)| 
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Proof of Eq. (HI: 



(G(''y-(rm) 



q(0) _ p(0) 



(G(o))^G(<') 



_ (^r(^-i))^G(^"') + --- + (r^^'j^G'"' 

+ (•plA'-ljj^pliV-l) ^ ^ (p(0)-)rp(0) 

= cg'^^^'^^g'""^' 



^-2 
^0 



g('v-i))^ rr(^-i) _^ j(A'-i)r(^-2) + . . . + j(^-i) . . . j(i)r(o)" 

.. + (r(o)) 



N-l 



«=0 

(G(^-i))^MG(^-i)-(G(^-i))^n-'£'[(rW)V")] 

«=o 



(57) 



For all j such that F-"^ ^ 0, the corresponding yj"' = and for all / such that F^-"'^ = 0, the coiTesponding 7;- ^ for 
all n. Thus T^") is always orthogonal to F'"'; therefore l^'^-^ [(r(''))^F(")] = 0. Thus, beginning with Eq. ^ and 
using Eq. ( flTl l, 

= (G(^-iy[MG(^-i)-n 
= a(g(^-^))V^-') 



= XF' 



(58) 



Here we have used Eqs. © and O to write = A2(G(^-^))^G(^-1). 
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